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Abstract. In this paper we consider the global existence of weak 
solutions to a class of Quantum Hydrodynamics (QHD) systems 
with initial data, arbitrarily large in the energy norm. These type 
of models, initially proposed by Madelung [ID], have been exten- 
sively used in Physics to investigate Supefluidity and Supercon- 
ductivity phenomena [19], [35] and more recently in the modeling 
of semiconductor devices [20] . Our approach is based on var- 
ious tools, namely the wave functions polar decomposition, the 
construction of approximate solution via a fractional steps method 
which iterates a Schrodinger Madelung picture with a suitable wave 
function updating mechanism. Therefore several a priori bounds 
of energy, dispersive and local smoothing type allow us to prove 
the compactness of the approximating sequences. No uniqueness 
result is provided. 



1. INTRODUCTION 

In this paper we study the Cauchy problem for the Quantum Hy- 
drodynamics (QHD) system: 

(1) , 

dtp + div J = 

d^j + div + VP(p) + pVV + /(Vp, J, VVp) = f pV 

-Ay = p-C(x), 

with initial data 

(2) p(0) = po, J(0) = Jo. 

We are interested to study the global existence in the class of finite 
energy initial data without higher regularity hypotheses or smallness 
assumptions. 

The analysis of uniqueness of weak solutions in some restricted classes 
will be done in a forthcoming paper. 

We will discuss in particular the case /(i/p, J, V^/p) = J, however we 
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are able to treat a more general collision term, as it will be clarified in 
the following Remark O Indeed, we show that it is possible to consider 
a general collision term of the form f = aJ + pVg, where a > and g 
is a nonlinear operator of ^/p, J, V^p, satisfying certain Caratheodory- 
type conditions (see Remark [5] for more precise explanations). 
There is a formal analogy between ([T]) and the classical fluid mechan- 
ics system, in particular when h = 0, the system ([T]) formally coincides 
with the (nonhomogeneous) Euler-Poisson incompressible fluid system. 
The theoretical description of microphysical systems is generally based 
on the wave mechanics of Schrodinger , the matrix mechanics of Heisen- 
berg or the path- integral mechanics of Feynman (see [19]). Another ap- 
proach to Quantum mechanics was taken by Madelung and de Broglie 
(see [IQ]), in particular the hydrodynamic theory of quantum mechan- 
ics has been later extended by de Broglie (the idea of "double solution" ) 
and used as a scheme for quasicausal interpretation of microphysical 
systems. There is an extensive literature (see for example [I2], [35] . 
[32] . [3T] . [M] . [30], and references therein) where superfluidity phe- 
nomena are described by means of quantum hydrodynamic systems. 
Furthermore, the Quantum Hydrodynamics system is well known in 
literature since it has been used in modeling semiconductor devices 
at nanometric scales (see [20]). The hydrodynamical formulation for 
quantum mechanics is quite useful with respect to other descriptions 
for semiconductor devices, such as those based on Wigner-Poisson or 
Schrodinger-Poisson, since kinetic or Schrodinger equations are com- 
putationally very expensive. For a derivation of the QHD system we 
refer to [2], [B], [2B], M, IS], US], M- 

The unknowns p, J represent the charge and the current densities re- 
spectively, -P(p) the classical pressure which we assume to satisfy P{p) = 
l^p^ (here and throughout the paper we assume I < p < 5). The 
function V is the self-consistent electric potential, given by the Poisson 
equation, the function C{x) represents the density of the background 
positively charged ions. In the paper we will treat extensively the case 
C{x) = 0, but all the results can be extended to more general cases. 
For instance we can assume C G W^^'-^(]R^) fl 14^^'^(M^). 



tential, or as a quantum correction to the pressure, indeed with some 
regularity assumptions we can write the dispersive term in different 
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can be interpreted as the quantum Bohm po- 
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ways: 

= —AVp-h'^div{Vy/p®Vy/p). 

There is a formal equivalence between the system ([1]) and the following 
nonlinear Schrodinger-Poisson system: 

. .^ / ihdtij + f = + vij + vij 

^ ^ \ -AV=\^\^ 

where V = log ^= j , in particular the hydrodynamic system ([T]) can 

be obtained by defining p = \ip\'^, J = him ('?/'VV') and by computing 
the related balance laws. 

This problem has to face a serious mathematical difficulty connected 
with the need to solve (jlj) with the ill-posed potential V. Presently 
there are no mathematical results concerning the solutions to (jlj), ex- 
cept small perturbations around constant plane waves or local existence 
results under various severe restrictions (see [27], [36]). In the paper 
by Li and the second author [37] there is a global existence result for 
the system ([1]), regarding small perturbations in higher Sobolev norms 
of subsonic stationary solutions, with periodic boundary conditions. 
A possible way to circumvent this type of difficulty could be to develop 
a theory regarding wave functions taking values on Riemann manifolds 
but we will not pursue this direction in this paper (see also [TU], [46j). 
Another nontrivial problem concerning the derivation of solutions to 
([1]) starting from the solutions to (jl]), regards the reconstruction of the 
initial datum ip^O) in terms of the observables p(0), J(0). Actually this 
is a case of a more general important problem in physics, pointed out 
by Weigert in [53j- He named it the Pauli problem (since this question 
originated from a footnote in Pauli's article in Handbuch der Physik, 
see jlS]), and it regards the possibility of reconstructing a pure quan- 
tum state, just by knowing a finite set of measurements of the state (in 
our case, the mass and current densities). Here the possible existence 
of nodal regions, or vacuum in fluid terms, namely where p = 0, forbids 
in general this reconstruction in a classical way, and in any case some 
additional requirements (quantization rules like the Bohr-Sommerfeld 
rule) would be necessary. In any case, various authors (see [53] and 
references therein) showed that knowledge of only position and mo- 
mentum distribution does not specify any single state. 
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The opposite direction, namely the derivation of solutions of (jl]) start- 
ing from solutions of ([T]) also can face severe mathematical difficulties 
in various points. In particular if we prescribe ^^(O), we can define p(0) 
and J(0), however from the evolution of the quantities p(t) and Jit), 
we cannot reconstruct the wave function ip{t). Furthermore, from the 
moment equation in ([1]) we cannot derive the quantum eikonal equation 

(5) d,S + l|V^r + hip) + V + S= 

which is the key element to reconstruct a solution of (jlj) via the WKB 
ansatz for the wave function, ip = ^e^^l^. 

Similar difficulties arise when approaching with Wigner functions, which 
has been recently quite popular to deduce quantum fluid systems in 
a kinetic way (see [Ml)- Even in the case we know the initial data 
p(0), J(0) to be originated from a wave function '0(0), it is very diffi- 
cult to show that the solutions p(i), J it) coincide for all times with the 
first and second momenta of the Wigner function obtained by solving 
the Wigner quantum transport equation. Moreover in our case there 
is also the difficulty due to the non-classical potential V . 
A related question has been investigated by Bourgain, Brezis, Mironescu 
in several papers (see [6] and references therein) regarding the lifting 
problem for harmonic maps. 

Another formal approach is provided by the following transport equa- 
tion 

(6) pdtS + jv^ + phip) ^pv^ps = yVp^Vp 

which is obtained by multiplying equation ([5]) by p. Unfortunately the 
lack of regularity of the solutions does not allow to apply even the 
more recent advances of the theory of transport equations [16], [1] or 
the somehow related approach developed by Teufel and Tumulka [51] 
in the study of trajectories of Bohmian mechanics. 
A related problem arises in the study of Nelson stochastic mechan- 
ics (see for instance Nelson [12], Guerra and Morato [24j), where the 
mathematical theory is based on the analysis of the velocity fields (see 
Carlen [7]). The application of this approach to our problem presents 
the same level of difficulty of the previously mentioned methods from 
the transport theory. 

A natural framework to study the existence of the weak solutions to 
is given by the space of finite energy states. Here the energy associated 
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to the system ([1]) is given by 

(7) E{t):= [ ^\V^)\' + l\A{t)\' + f{p{t)) + hvV{t)\Mx, 

where A := J/ y/p, and /(p) = ■:^p^ ■ The function /(p) denotes the 
internal energy, which is related to the pressure through the identity 
P{p)=pf{p)-f{p). 

Therefore our initial data are required to satisfy 

(8) Eo := / %\VVp'^? + ^|Ao|' + /(po) + ^|Vl^o|'dx < oo, 

or equivalently (if we have 1 < p < 5), 

(9) v^G//^(M3) Ao:= Jo/v/p^GL2(M3). 

Definition 1. We say the pair (p, J) is a weak solution of the Cauchy 
problem ([1]), ([2]) in [0, T) x with Cauchy data (po, Jo) e ^^^(^3)^ 
there exist locally integrable functions A, such that ^/p G Lf^cil'^^ ^ 
A G LL([0,T);LL(M3)) and by defining p := {^f , J := ^A, one 
has 

• for any test function rj G C^([0,T) x M.^) we have 

(10) / / pdtri + J ■ Vrjdxdt + [ po??(0)dx = 0; 



• for any test function C G C^([0, T) x M^; RS) 

(11) / / J-9iC + A®A: VC + P(p)divC-pVr-C- J-C 

Jo Jr^ 

+ ffVy/p ® V^/p : VC - -rpA div Cdxdt + Jq ■ C(0)dx = 0; 

• (generalized irrotationality condition) for almost every t G (0, T), 

(12) VAJ = 2Vv/pAA 
/ioWs m t/ie sense of distributions. 

Remark 2. Suppose we are in the smooth case, so that we can factorize 
J = pu, for some current velocity field u, then the last condition ( fT2l) 
simply means pV Am = 0, the current velocity u is irrotational in pdx. 
This is why we will call it generalized irrotationality condition. 

Definition 3. We say that the weak solution (p, J) to the Cauchy prob- 
lem (Q, is a finite energy weak solution (FEWS) in [0,T) x M^, if 
in addition for almost every t G [0,T), the energy ^ is finite. 
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In the sequel we restrict our attention only to FEWS with po ^ 
L^(M^). The hydrodynamic structure of the system ([T]), ([2]) should not 
lead to conclude that the solutions behave like classical fluids. Indeed 
the connection with Schrodinger equations suggests that po, Jq should 
in any case be seen as momenta related to some wave function ipQ. The 
main result of this paper is the existence of FEWS by assuming the 
initial data po, Jo are momenta of some wave function ipQ G iJ^(]R^). 

Theorem 4 (Main Theorem). Let ipo € H^{M.^) and let us define 

Po ■= l^oP, Jo ■= Mm(V'oV?/'o)- 

Then, for each < T < oo, there exists a finite energy weak solution 
to the QHD system ([1]) in [0,T) x M^, with initial data {po, Jo) defined 
as above. 

Remark 5. As we said above, we can extend the result of the previous 
main Theorem \^to a more general type of collision term. Indeed it 
is possible to consider nonlinear terms of the type f{^/PiJi'^^/p) = 
aJ + pVg(t, X, ^/p, A, V^), where a > and g satisfies the following 
Caratheodory type conditions 

• for all {^/p, A) G [0,oo)xR^, the function {t,x) t-^ g{t, x, y/p, A,V y/p) 
is Lebesgue measurable; 

• for almost all {t,x) G [0, oo) x M^, the function {u,v,w) i— >■ 
g(t,x,u,v,w) is continuous for {u,v,w) G [0, cxd) x x M^; 

• there exists C G L°°{\Q, oo) x R^) such that 

(13) \g{t,x,u,v,w)\ < C(t,x)(l + \u\'^ +\v\^l^ + \w\^'^). 

The case a = requires a slight modification of the method presented 
in the Section^ as we are going to explain below. The case a > can 
be treated with minor modifications of the finite difference scheme in 
the Section\^ defined for f{y/p, J, V^p) = J. See the Remark\2^ below 
in the sectionlB 

We remark that the condition (|T3l) is needed in order to ensure g G 
L°°([0,oo);Lp(R3) + L°°(R3)) as long as ^ e L°°([0, T); //^(R^)), 
A G L°°([0,T);L2(R3)). Note that with this definition of f we need 
to change the definition given before for the weak solutions. Indeed the 
balance law for the current density becomes 




J ■ + A ® A : VC + P(p) divC - pVV -C-J-C 



+ g{t,x, v/p,A)(Vp-C + pdivC) 

fi^ r 

+ h^Vy/p ® V Vp : VC - — pA div Cdxdt + Jq ■ C(0)da; = 0. 

4 Jr3 
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In the Sections and we will explain how to modify our methods to 
include this case (see Remarks WE and\EW- 

At the end of this Section we would hke to remark that the system 
([1]) is closely related to the Quantum Drift-Diffusion (QDD) equation 

(14) dtp + div (PV (y ^ - V - Mp)) ) , 

where V is again the electrostatic potential and the enthalpy h{p) is 
such that ph'{p) = P'{p). Indeed if we scale the collision term J in ([1]) 
with a relaxation time e and we write 
(15) 

{dtp' + div J' = 
dtJ' + div + VP(p^) + p'W + \J' = f p^V (^) 

-AV = p'' -C{x), 

then in the formal limit e — we get the equation f|T^ . A partial 
result on this limit was given by Jiingel, Li and Matsumura in ^6] . 
while in [22] it is shown the global existence of non-negative variational 
solutions to ( fT^ without the electrostatic potential and the pressure 
term. In a forthcoming paper we deal with the relaxation limit from 
solutions to (fTS!) to solutions of (fT^ . 



2. Preliminaries and Notations 

2.1. Notations. For convenience of the reader we will set some nota- 
tions which will be used in the sequel. 

If X, Y are two quantities (typically non-negative), we use X <Y to 
denote X < CY, for some absolute constant C > 0. 
We will use the standard Lebesgue norms for complex-valued measur- 
able functions / : M°' ^ C 

1 



If we replace C by a Banach space X, we will adopt the notation 

_ i/p 

\lp{R'';X) ■= ( / \\f{x)\\xdx 

to denote the norm of / : M'^ — > X. In particular, if X is a Lebesgue 
space L^(M'^), and d = 1, we will shorten the notation by writing 

L?Ls(/xR") :=( [ \\fmU^^)dt) =( [{[ \f{t,x)\^dxy/^dt^ 
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to denote the mixed Lebesgue norm of /:/—>■ L'' 
will write LlLl{I x M") := L'?(/; L''{R'^)). 

For s G M we will define the Sobolev space H'{W) := {l-A)-'/^L'^{W). 

Definition 6 (see p]). We say that {q,r) is an admissible pair of 

exponents if 2 < q < oo, 2 < r < Q, and 

13/11 

Now we will introduce the Strichartz norms. For more details, we 
refer the reader to [50], [11]. Let J x be a space-time slab, we define 
the Strichartz norm 5'°(/ x M^) 

X 1/2 

I'^IUoaxMS) •= sup I ^ l|-PAf^llL9LJ{7xI 
V N 

where the sup is taken over all the admissible pairs {q,r). Here 
denotes the Paley-Littlewood projection operator, with the sum taken 
over diadic numbers of the form N = 2\ j G Z. For any A; > 1, we can 
define 

Note that, from Paley-Littlewood inequality we have 

\\u\\LlLr(IxR3) < II (^ \PNu\'^y^^\\LlLriIxR^) < Yl H^^^Hi^ 
N \ N 

and hence for each admissible pair of exponents, one has 

II'"IIl?LJ(/xR") ^ ll^lls-oC/xRS)- 

Lemma 7 (p^). 

ll'"llL4Lg°(/xIR3) < ||'"|l5i(/xM3)- 

2.2. Schrodinger equations. In this paragraph we recall some im- 
portant results concerning the nonlinear Schrodinger equations that 
will be used throughout the paper. First of all, we recall a global well- 
posedness theorem for the nonlinear Schrodinger-Poisson system (see 
[8] and the references therein): 

Theorem 8. Let us consider the following nonlinear Schrodinger-Poisson 
system 

1. 7^ f ihdtiJ + f = |^|P- V + Vi) 



1/2 
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with the initial datum 

(18) i/j{0)=^oeH\R'). 

There exists a unique globally defined strong solution ip G C(R; H^(M.^)), 
which depends continuously on the initial data. Furthermore, the total 
energy 

(19) E{t):= [ ^\V^P{t,x)\' + ^m,x)r' + V{t,x)m,x)\'dx 

Jm^ 2 p+1 

is conserved, namely 

(20) E{t) = Eo, for all t G M. 

Moreover, the dispersive nature of Schrodinger equation provides also 
some further integrabihty and regularity properties of the solutions. 
The first result in this direction regards some space-time integrabihty 
properties, the most important of them are the well known Strichartz 
estimates for the Schrodinger equation in R^+^. We refer to the classical 
paper of Ginibre and Velo [23], the paper of Keel and Tao [33j, the 
books of Cazenave [8] and Tao [50] and the references therein. Let 
ip is the unique (strong) solution of the Cauchy problem for the free 
Schrodinger equation 

r idt^ + = 

I m = ^0, 

in the following we shalll denote by f/(-) the free Schrodinger group, 
defined by the indentity U{t)ipQ = ipit). 

The next result is taken from [33] and the estimates are obtained in a 
very general setting. 

Theorem 9 (Keel, Tao [35). Let [q, r), (g, f) be two arbitrary admissi- 
ble pairs of exponents, and let U{-) be the free Schrodinger group. Then 
we have 



(21) \\Uit)f\\,.,r < ||/|U.(M3) 

(22) II / Uit-s)Fis)dsh.Lr<\\F\\^,^„ 

J s<t * 

(23) II j U{t - s)Fis)ds\\r^.Lr < II^^ILf ^j'- 



By using the previous Theorem, we can state the following Lemma 
about some further integrabihty properties for solutions of the Schrodinger 
equations. 
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Lemma 10 ([llj). Let I be a compact interval, and letu:IxM.^^C 
be a Schwartz solution to the Schrddinger equation 

idtu + Au = Fi + ... + Fm, 

for some Schwartz functions Fi, . . . , Fm- Then we have 

where (gi, ri), . . . , {qm, tm) are arbitrary admissible pairs of exponents. 

Furthermore, the solutions to the Schrodinger equation enjoy some 
nice local smoothing properties; there are various results in the litera- 
ture regarding this property, here we recall a theorem, due to Constatin 
and Saut [9j which actually covers a more general setting. In particular 
it can been shown the solutions of the Schrodinger equation are locally 
more regular than their initial data. Close results have been proved by 
P. Sjolin [38] and L. Vega p]. 



(24) 



Theorem 11 (Constantin, Saut [9]). Letu solves the free Schrddinger 
equation 

idtu + Am = 
m(0) = Uq. 

Let X e C^(Ki+3) of the form 

Xit,x) = Xo{t)xi{xi)X2ix2)x{x3), 

with Xj £ Cq^ (M) ■ Then we have 

f x\t,x)\{I - Ay/\{t,x)\Mxdt < C^WuoWl'^us. 

In particular, if Uq G L^(]R^), one has for all T > 

ueL'{[0,T];Hli\R')). 

We have a similar result also for the nonhomogeneous case: 

Theorem 12 (Constantin, Saut [9J). Let u be the solution of 

idtU + A?i = F 
m(0) = Mo G L2(M3)^ 

where F e L^{[0,T]; L'^{R^)). Then %t follows 

ueL\[0,T]-Hli\R')). 

Moreover, let x £ C^(M^"^^) be of the form 

X{t,x) = Xo{t)xi{xi)X2{x2)X3{x3) 



(25) 
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with Xj £ ^o°(-'^)' suppxo C [0,T]. Then the following local smoothing 
estimate holds 

(26) ([ x\t,x)\{I-A)'/\it,x)\Mxdty 

< C (||Mo||l2(ir3) + ||-^||l1l2([o,t]xR3)) 

This results imply that, the free Schrodinger group U{-) fulfills the 
following inequalities 

(27) \\U{-) 

(28) II / f/(t - s)F(s)ds||^2(jQ^^j.^v2(^3^) < ||F||li([o,t];L2{r3)). 

^ 

2.3. Compactness tools. Here we recall some compactness theorems 
in function spaces, which will be relevant in the Section [6] to prove the 
convergence of the approximate solutions. Let us recall in particular 
a compactness result due to Rakotoson, Temam [H] in the spirit of 
classical results of Aubin, Lions and Simon, see [17], [3], [39] . 

Theorem 13 (Rakotoson, Temam 01]). Let {V, \\ ■ ||y), {H; \\ ■ \\h) be 
two separable Hilbert spaces. Assume that V (Z H with a compact and 
dense embedding. Consider a sequence {m^}, converging weakly to a 
function u in L^([0,T]; V), T < oo. Then -u^ converges strongly to u 
in -^v^([0, T]; if ), if and only if 

(1) M^(t) converges to u{t) weakly in H for a.e. t; 

(2) lim|i5|_,o,i?c[o,T]Sup^>o/ij h%t)\\jjdt = 0. 

3. Polar Decomposition 

In this section we will explain how to decompose an arbitrary wave 
function ip into its amplitude ^/p = \^jJ\ and its unitary factor 0, namely 
a function taking its values in the unitary circle of the complex plane, 
such that ijj = ^Jp(^. The idea is similar in the spirit, to that one 
used by Y. Brenier in |5], to find the measure preserving maps needed 
to write vector-valued functions as compositions of gradients of con- 
vex functions and measure preserving maps. Our case is much simpler 
than [5j and it can be studied directly in a simpler setting. Brenier's 
idea looks for projections of functions u onto the set of measure 
preserving maps S", contained in a given sphere of L^, i.e. one has to 
find s G S", which minimizes the distance ||ti — s 11^2, or equivalently, 
which maximizes (m, 5)^2. In our case we should maximize Re(M, 5)^2, 
within complex-valued functions with the constraint to take value in 
the unit ball of Actually in this case the maximizer turns out to 
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be trivially determined and the previous variational argument has only 
to be considered a motivation to the subsequent considerations. 
However in the case we assume u G H-^lM.^) (that is the relevant situa- 
tion in physics), then it would be useful to know the regularity of the 
maximizers in the light of a possible use of the considerations explained 
in the Remark [271 

Let us consider a wave function ip G L^(]R^) and define the set 
(29) 

P{ip) := {4> measurable | ||0||l°°(r3) < 1, = ip a. e. in M'^} , 

where ^ = IV'I- Of course if we consider G Pijp), then by the 
definition of the set P{ip) it is immediate that \(j)\ = 1 a.e. ^/p — dx in 
and G P{ip) is uniquely determined a.e. — da; in M^. 

Remark 14. Let B^i be the ball in with radius R, centered at the 
origin and let us define the set 

PrU') ■= {0 measurable \ 00(^3) < 1, ^/p(|) = ip a. e. in M^} . 

Then it is easy to see that G Pniip) if o-nd only if (p is a maximizer 
for the functional 

$i?[0] := Re / %lj(t)dx 

over the set 

Sr := {0 G L^Bn) \ UWl-^Br) < l} • 



The next lemma connects the structure of the bilinear term h'^KeCVip ® Vip) 
with V^/p ® V^/p and Moreover it shows this structure is H^- 

stable. 



Lemma 15. Let tp G H^{^), := 1-01, then there exists G 
such that Ip = yp0 a.e. in R^, y/p G H^(M.^), V^/p = Re(0VV'). If we 
set A := hlm{(pVip) , one has A G L^(M'^) and moreover the following 
identity holds 

(30) ffRe{dji^dki^) = ffdj^dk^ + A^-'U^*^). 
Furthermore, let tpn ^ ip strongly in H^{M.^), then it follows 

(31) Vv/p;^VVp, A„^A inL^W"), 
where A„ := /iIm(0„V'0n)- 

Proof. Let us consider a sequence {ipn} C C°°(M^), -0 — -0 in H^{M^^), 
define as before 

(32) 0„(x) = <^ IV-nWI -^"^^ ^ 

[ if '0(x) = 0. 
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Then, there exists (p G L°°(M3) such that 0„ ^ in L°°{R^) and 
Vijjn Vijj in L2(M3), hence 

Re (J^Vipn) Re (0V^) in L'^(M.^). 

Since by fl32|) . one has 

Re (0nVV^„) = V|?/'„| 

a.e. in M^, it follows 

V ^ Re (0VV^) in (M^) . 
Moreover, one has ^ in L^(]R^), therefore 

Vv/;^ = Re (0V^) , 
where is a unitary factor of ip. 

The identity fl30|) follows immediately fro the following 

n2Re(ajV^4V) = h'^Re{{(j)dji^)(^dk^)) = ffRe{(j)dj'^)Re(^dk^lj) 

Now we prove (13T1) . Let us take a sequence ipn ^ '4' strongly in 
and consider V ^/p^ = Re (0„VV'n), A„ := /ilm (0„V^/'„). As before, 
0„ ^0 in L-(R3), with a polar factor for_^; then - V^/p, 

Re (0^VV^„) ^ Re (^VV'), and V^/p = Re (^VV')- Moreover, A„ : = 
Mm [(pn^tpn) Mill (^V^/') =: A. To upgrade the weak convergence 
into the strong one, simply notice that by fl30|) one has 

^2||v^||2, = ff\\V^\\l, + \\A\\l, < liminf (/i^H Vy/p^Hi, + ||A„||i.) 

n—^oa 

= liminf n^ll V^„||i2 = n2||V^||i2. 

n— >oo 

□ 



Corollary 16. Let ip G H^W), then 

(33) A VV' = 2iVv/p A A. 

Proof. It suffices to note that we can write 

WV' A VV' = h{(j)Vi)) A (0V^), 

where G L°°(R^) is the polar factor of such that V^/p = Re(0VV'), 
A := Mm(0VV^), as in the Lemma [T5l By splitting into its real 
and imaginary part, we get the identity (!33|l . □ 
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Now we state a technical lemma which will be used in the next sections. 
It summarizes the results of this section we use later on and will be 
handy for the application to the fractional step method. 

Lemma 17. Let ip G H'^iM?), and let r, e > he two arbitrary (small) 
real numbers. Then there exists ip G H^{M.^) such that 

P = P 

A = (l-r)A + r„ 

where := Itp], y/p := \'ip\,-A '■= hlm{(f)Vijj), A := hlm{(f)V4') , 0, are 
polar factors for ip , ip , respectively, and 

Furthermore we have 

(34) VV^ = VV^ - i^^'^A + r,,,, 
where ||0*||j;^oo(ir3) < 1 and 

(35) WrsAL^R^) < C(^|| V^/'||l2(b;3) + e). 

Proof. Consider a sequence {ipn} C C°°(M^) converging to ip in if^(M^), 
and define 

^ / ^ f rH4 ifV^n(x) ^0 
0n(x):= I'^-WI ^ 

[ if ijjnix) = 

as polar factors for the wave functions ipn- Since ^pn ^ then 0„ is 
piecewise smooth, and fin '■= {a; G : |^/'n(a^)| > 0} is an open set, 
with smooth boundary. Therefore we can say there exists a function 
9n '■ fin [0, 27r), piecewise smooth in fin and 

0„(x) = e*^"*^^\ for any x e Qn- 



Moreover, by the previous Lemma, we have 0„ ^ in L°°(M ), where 
is a polar factor of tp, and A„ := hlm{(pnVipn) —>■ A := hlm{(p'Vip) in 
L^fM^). Thus there exists n eN such that 



11^ - ^n||/fi(K=^) + 11^ " ^n\\L2(R-i) < E. 

Now we can define 

(36) ;p:=e'^^-^^'-^n. 



Let us remark that (!36|) is a good definition for -0, since outside fl^ 
where On is not defined, we have = 0, and hence also ipn = 
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Furthermore, 



n 



+ (VV^„ - VV^) - z^e*(i-^)^"(A„ - A) 
n 

and thus we obtain (IMl) and ( |35l) . with r^^r given by 

r.,. = r ^^^^^^ j + (V^„ - V^) - z^e^(^-^)^"(A„ - A) 

Moreover 

(37) A = Mm(e^*(^-")^"VV^) = (1-t)A„ = (1 -r)A + (1 -r)(A„ - A) 

and clearly := (1 — r)(A„ — A) has L^-norm less than e by assump- 
tion. □ 



4. QHD WITHOUT COLLISIONS 

Now let us summarize some key points regarding the existence of 
weak solutions to the Quantum Hydrodynamic system, in the colli- 
sionless case. 

The balance equations can be written in the following way: 

{dtp + div J = 
dtJ + div (:^) + VP(p) + pVV = f pV (^) 
-AV = p, 

where P(p) = £^p(P+^)/^ 1 < p < 5. 

Definition 18. We say the pair (p, J) is a weak solution in [0, T) x 
M'^ to t/ie system (!38|) wzt/i Cauchy data (po, Jo) ^ -^Lc(^^) ^/ ^'^^ 
on/y z/ t/iere exzsi locally integrable functions ^/p, A such that y/p G 
LL([0,T);/7i,,(M3)),A e LL([0, T); LLM) an^^ z/ w;e de/^ne p : = 

{y/pY> J •= v^"^' ^^^^ 

• /or any test function G C(f ([0,T) x M^) we have 

(39) /" /" pdt^ + J ■ Vy^dxdt + f po¥5(0)dx = 

Jo Jrs Jrs 
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(41) 



and for any test function t] G C^([0, T) x M'^; M'^) we have 

(40) [ [ J-dtr] + A^A:Vr] + P{p)divr]-pW-rj 
Jo Jr^ 

r 

+ h'^Vy/p ® Vy/p : Vt] - —pA div r]dxdt + Jq ■ r]{0)dx = 0; 

• (generalized irrotationality condition) for almost every t G (0, T) 
V A J = 2Vv/p A A 

/ioWs m t/ie sense of distributions. 

We say that (p, J) zs a finite energy weak solution if and only if it is 
a weak solution and the energy E{t) is finite a.e. in [0,T). 

The next existence result roughly speaking shows how to get a weak 
solution to the system (l38l) out of a strong solution to the Schrodinger- 
Poisson system 

ihdtip + ^Aip = + Vi/j 

-AV= IV'P 

The quadratic nonlinearities in fl38l) are originated by a term of the 
form Re {Vip ® ^ip) since formally 

h^Re (V^ ® Vip) = h'^Vy/p ® V^/p + i^-L. 

However this identity can be justified in the nodal region {p = 0} only 
by means of the polar factorization discussed in the previous section. 
Indeed, we stress that in the previous identity the right hand side is 
written in terms of p and J which do exist in the whole space M^, via 
the Madelung transform. However, the term should be interpreted 
as A® A, where A is the Radon- Nykodim derivative of Jdx with respect 
to y/pdx. Unfortunately the Madelung transformations are unable to 
define A on the whole M^, hence one needs to use the polar factorization 
discussed in the previous section to define A in the whole M^. 
Furthermore, the study of the existence of weak solutions of fl38l) is 
done with Cauchy data of the form {po, Jq) = {\ipo\'^,hlm{xljQ'V4'o)), 
for some tpo G if^(M^;C). These special initial data yield to consider 
the Cauchy problem for the QHD system only for solutions compatible 
with a wave mechanics point of view. 



Proposition 19. Let < T < oo, let ipo e if (M"^) and define the 
initial data for ( |38|) . (po, -^o) '■= (IV'oP; Klm^ipoVipo)) ■ Then there exists 
a finite energy weak solution (p, J) to the Cauchy problem ( 1381) in the 
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space-time slab [0, T) x M^. Furthermore the energy E{t) defined in (JTj) 
is conserved for all times t G [0, T). 

The idea behind the proof of this Proposition is the following. Let 
us consider the Schrodinger-Poisson system (14T]) with initial datum 
= ipQ. It is well known that it is globally well-posed for initial 
data in H^(R^) (see [8j), and the solution satisfies G C°(M; H\R^)). 
Then it makes sense to define for each time t G [0,T) the quantities 
p{t) := |'?/'(t)p, J(t) := hlm{ip{t)Vip(t)) and we can see that [p, J) is a 
finite energy weak solution of (!38l) : indeed for the solution ip of (l4TI) 
we have 

in the sense of distributions. Thus formally we get the following iden- 
tities for (p, J) 

dtp = —hlm(^pAijj) = — div (^^^^{ipVip)) , 



dtJ = him (vi^ |^-yA^+ -^(1^1^-1 + V)'^^^ 
+ him (ip (yAV^ - -^Vd^r^ + V)tP - ^M"'' + V)VtP^^ 

= ^VAl^jjf - h^ div (Re(V^ ® V^)) - ^^Vd^/'l^'+i) - I^I^VK 

Thanks to the Lemma [15] we can write 

^^Re (VV^ ® V?A) = h^V^/p ® + A ® A 

with y/p, A as in statement of the Lemma [151 ^-iid clearly they also 
satisfy 

J= V^A. 

Hence formally the following identity holds: 

h"^ 

(42) 9tJ+div(A®A)+VP(p)+pVy = — AVp-^Miv(Vv/p®Vv/p). 

Of course these calculations are just formal, since ip doesn't have the 
necessary regularity to implement them. 
Furthermore, it is well known that the energy 
(43) 

m--= [ ^l^i'l'+^li'r'dx+l [ \^p(t,y)\'^^Mt,x)\'dxdy 
V 2 p+1 2 JmxRs \x-y\ 
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is conserved for the Schrodinger-Poisson system ( HTi) . Thus, it only 
remains to note that by Lemma [13 the energy in is equal that one 
in ©. 

Proof. Let us consider the Schrodinger-Poisson system (HTi). From stan- 
dard theory about nonlinear Schrodinger equations it is well known 
that pij) is globally well-posed for initial data in the space of en- 
ergy: V(0) = ^0 £ if"'^(M^). Now, let us take a sequence of mollifiers 
{Xe} converging to the Dirac mass and define ip'^ := Xe * where 
i(j e C{R;H^{M.^)) is the solution to dH]). Then tp^ e C°°{R^+^) and 
moreover 

(44) ihdti^' + = xe^ (iv-r-v + vijj) , 

where V is the classical Hartree potential. Therefore differentiating 
lip^l"^ with respect to time we get 

dtin^ = 2Reiifdtr) = 2Re (^^(yA^^ - ^Xe * (1^1^"^ + Vt/j)) 

= -hdiv {lm{ifVip')) + ^Im {i^Xe * (I^T'V + Vip)) . 
If we differentiate with respect to time hlmlip^'Vip'^) we get 



dt (Mm(^^VV'^)) = him 



+ him 

h^ 



2 Re (-V^Z-'A^ + ^AV^')+Re {xe * (|^r"^^)V^" - i^Xe * V 

+ Re (x ^ {V^W - * V(\/V')V^) =:A + B + C 

Let us discuss these three terms separately. For the first term A it is 
immediate that 

^2 fc2 

yRe {-Vij'Ai:^ + ^j^AVi^') = h^ div (Re(V^^ ® Vi)'))+—AV\i)'\^. 

Let us discuss the second and the third terms, B and C. We can recast 
B in the following way 

Xe ^ - T'xe * vdv^r- V) 

= \ij'\p-^^vij' - ^xe * (V^v|^|p-^ + + 
= -i^Xe* ii^'^W-') +R{ + Rl = -\r?^\rr' + RI + RI + Rl 
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where 

Ri := {xe ^ m"-'^) - m^'W) 
Ri := {\rr'vr -xe^ (i^r v^)) 
Ri := r {ry\rr' - ^ (v^vi^r-^)) . 

Now it only remains to analyze how the remainder terms i?^ go to 

zero. We will prove it in the next Lemma. Furthermore, the following 
identity holds 

' ' ' ' p+ 1 ' ' 

For the third term C after similar computations we get 

C=mV{xe^V) + Rl, 

where the remainder term R^ will be analyzed in the next Lemma. 
Hence we can conclude that hlm{ilj'^'Vip'^) satisfies 

p+i 

where R^ := i?f + + i?| + RI- Now, as e 0, we get 

- - 

dt (Mm(V^VV^)) = -h^ div (Re(V^ ® V^)) + —^V\ip\'^ 

This identity is equivalent to (H2|) since by the Lemma [15] we have 

^^Re {Vi) ® Vip) = f?^4p ® V^/p + A A, 
where A is defined as in the Lemma [15] and y^A = J, where J := 
Mm(?/'V'?/'). 

Finally, let us note that, by the definition of J, we have 

V A J = himiyi) A VV") 

then by the Corollary [16], we get V A J = 2V^ A A. □ 



Lemma 20. Let < T < oo, then 

||-Rj||Li^([o,T]xR3) ^0, as e ^ 

J = 1,2,3,4. 
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Proof. It is a direct consequence of Strichartz estimates for the solution 
of the Schrodinger-Poisson system fHTl) 

(45) II^IIl«h/i.'-([o,t]xr3) < C{Eo, ||V'o||l2(r3),T), 

where (g, r) is a pair of admissible exponents and Eq is the initial energy 
for fHTl) . The first error can be controlled in the following way 

\\RI\\lIJ10,T]xR3) 

1^^^+'' L," ([0,T]xR3) Lf(f-i^LP+i([0,r]xR3) 

' / rP+l 

The remaining error terms can be computed in the same way. We re- 
mark that by using the Strichartz estimates it folUows ilj'V\tp\^~'^ lies in 

4(p+l) 

L'^'-'^ LP+\[0,T]xR^). □ 



5. THE FRACTIONAL STEP: DEFINITIONS AND CONSISTENCY 

In this section we make use of the results of the previous sections to 
construct a sequence of approximate solutions of the QHD system 

{dtp + div J = 
dj + div (^^) + VP(p) + pVV + J= f pV 
-AV = p ' 
with Cauchy data 

(47) (p(0),J(0)) = (po,Jo). 

Definition 21. Let < T < oo. We say {(p^, J^)} is a sequence of 
approximate solutions to the system ([T]) in [0,T) xM^, with initial data 
{po, Jo) £ Lloci^^)' if there exist locally integrable functions y/p^,A'^, 
such that ^ e LL([0, T); //^(M^)), G LL([0, T); LL(ffi=')) and 
we define p"^ := (-^/p^)^, J"^ := i/p^A^, t/ien 

• /or any ies^ function rj G C^([0, T) x R^) one has 

(48) /" f p''dtT] + ■Vr]dxdt+ [ por/(0)dx = o(l) 
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Jo 




05 T — > 0; 

• (generalized irrotationality condition) for almost every t e (0, T) 



Our fractional step method is based on the following simple idea. 
We split the evolution of our problem into two separate steps. Let us 
fix a (small) parameter r > then in the former step we solve a non- 
coUisional QHD problem, while in the latter one we solve the coUisional 
problem without QHD, and at this point we can start again with the 
non-collisional QHD problem, the main difficulty being the updating of 
the initial data at each time step. Indeed, as remarked in the previous 
section we are able to solve the non-collisional QHD only in the case of 
Cauchy data compatible with the Schrodinger picture. This restriction 
impose to reconstruct a wave function at each time step. 
Now we explain how to set up the fractional step procedure which gen- 
erates the approximate solutions. We first remark that, as in the pre- 
vious section, this method requires special type of initial data {po, Jq), 
namely we assume there exists ipo G iJ^(M'^), such that the hydrody- 
namic initial data are given via the Madelung transforms 



This assumption is physically relevant since it implies the compati- 
bility of our solutions to the QHD problem with the wave mechanics 
approach. The iteration procedure can be defined in this following way. 
First of all, we take r > 0, which will be the time mesh unit; therefore 
we define the approximate solutions in each time interval [kr, (A;-|- 1)t), 
for any integer A; > 0. 

At the first step, /c = 0, we solve the Cauchy problem for the Schrodinger- 
Poisson system 



we have 



V A = 2Vv^ A A^. 



Jo = /iIm('0oV'0o)- 



(50) 



V'"(0) = V'o 
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by looking for the restriction in [0, r) of the unique strong solution ip G 
C{^-H^{R^)) (see 0). Let us define := := Mm(^V^/'^). 

Then, as shown in the previous section, (p"^, J^) is a weak solution to 
the non-coUisional QHD system. Let us assume that we know ip'^ in 
the space-time slab [{k — l)r, kr) x M^, we want to set up a recursive 
method, hence we have to show how to define ip'^,p'^, J'^ in the strip 
[kT,{k + l)r). 

In order to take into account the presence of the collisional term f = J 
we update int = kr, namely we define ilj'^{kT+). The construction of 
■ip'^{kT+) will be done by means of the polar decomposition described 
in the Section [3l 

Let us apply the Lemma [T71 with = ip'^^kr—), e = r2~'^||'?/'o||Hi(R3), 
then we can define 

(51) ip^{kT+) = ^, 

by using the wave function ip defined in the Lemma [T71 Therefore we 
have 

(52) p^ikT+) = p^ikr-) 

(53) A^ikT+) = {l-T)A^{kT-) + Rk, 

where ||-Rfc||L2(ig3) < r2"*=||V^o||//i(K=') 

(54) V^^{kT+) = V^^{kT-) - i^(P*A^{kT~) + rk,r, 
with < 1 and 

\\rkA\L^ < C(r|| VV^^(A;r-)|| + T2-%M\mm) < . 

We then solve the Schrodinger-Poisson system with initial data ipiO) = 
ip'^{kT+). We define ip'^ in the time strip [kr, (fc + l)r) as the restriction 
of the Schrodinger-Poisson solution just found in [0,r), furthermore, 
we define := \ip'^\'^,J'^ := hlm{ip'^Vip^) as the solution of the non- 
coUisional QHD ll38l). 

With this procedure we can go on every time strip and then construct 
an approximate solution (p"^, J'^,V'^) of the QHD system. 

Remark 22. To cover the general case where f{^: J, ^^/p) = olJ + 
pVg{t,x, ^yp. A, Vy^), a >0, g satisfies the conditions in the Remark 
and one has a non-zero doping profile, we solve the following Cauchy 



FINITE ENERGY WEAK SOLUTIONS TO QHD 



23 



problem 



(55) 



ihdtip + ^Aip = \ip\P-^ilj + V^ + g^ 
-A\/ = \i)\^-C{x) 
^(0) = ^0 



The global well-posedness follows from the result in [8]. 

Theorem 23 (Consistency of the approximate solutions). Let us con- 
sider the sequence of approximate solutions {{p^ , J^)}t>o constructed 
via the fractional step method, and assume there exists < T < oo, 
VP e LL([0, T); HUR^)) and A G LL([0, T); LUR')), such that 



(56) 
(57) 



p^^^ znL2([0,T);ifL(R3)) 
A^^A znL2([0,T);LL(M')). 



Then the limit function {p,J), where as before J = y/pA, is a weak 
solution of the QHD system, with Cauchy data (po)<^)■ 

Proof. Let us plug the approximate solutions (p"^, J'^) in the weak for- 
mulation and let C e C(f ([0,T) X ]R3), then we have 

/'OO /• 

/ / j-.9iC + A^®A-: VC + P(p^)divC-p^Vl^^-C-^"-C 
+Wv^(g) VVp^: VC- — P^AdivCdxdt+ / Jq ■ C(0)dx 

;"(fc+l)r 

= ^ / / J" ■ + A" ® A" : vc + P(p") divc - p^vv^ -c-r -c 



fc=0 



: VC - — p^ A div Cdxdt + / Jq ■ C(0)da; 



°o r p{k+l)i 
IkT 



fc=0 



J"-Cda;dt+ / J"((A;+l)r-)-C((A;+l)r) 



+ / Jo ■ C(0)dx 



~r{kT) ■ C(A;r)dx 
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V/ / -r.cdxdt + J2 {r{kT-)-.r{kT))-akT)d 

p OO „ 

/ -J" ■ (dxdt + V r / J"(A;r-) ■ C(A;r)da; 

oo „ 

+ V / (l-r)v/;^^i?,-C(fcr)da; 



E 



°o /•(A;+1)t 

E 



k=0 



kr 



as r ^ 0. 



J"((A; + l)r-) ■ C((A; + l)r) - J"(t) ■ C(t)dxdt 



oo „ 

+ (1 - r) V / v^(A;r)i?, ■ C(A;r)dx 

fc=0 

= o(l) + 0(r) 



□ 



6. A PRIORI ESTIMATES AND CONVERGENCE 

In this section we obtain various a priori estimates necessary to show 
the compactness of the sequence of approximate solutions (p^, J^) in 
some appropriate function spaces. As we stated in the Theorem |23l we 
wish to prove the strong convergence of {\^} in L^^^([0,T); Hl^^{R^)) 
and that one of {A^} in Ll^{[0,T)] Ll^{R^)). To achieve this goal we 
use a compactness result in the class of the Aubin-Lions's type lemma, 
due to Rakotoson-Temam [H] (see the Section 2). The plan of this 
section is the following, first of all we get a discrete version of the (dissi- 
pative) energy inequality for the system ([T]), later we use the Strichartz 
estimates for Vip'^ by means of the formula ( l6Ti) below. Consequently 
via the Strichartz estimates and by using the local smoothing results 
of Theorems [11], [12], we deduce some further regularity properties of 
the sequence {'Vip'^}. This fact will be stated in the next Proposition 
[331 In this way it is possible to get the regularity properties of the 
sequence {Vip'^} which are needed to apply the Theorem [13] and hence 
to get the convergence for {a/p^ and {A"^}. 

Let us begin with the energy inequality. First of all note that if we 
have a sufficiently regular solution of the QHD system, one has 

(58) E{t) = - [ [ \A\Mxdt' + Eo, 

Jo Jm^ 
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where the energy is defined as in (^^. Now we would hke to find a 
discrete version of the energy dissipation for the approximate solutions. 

Proposition 24 (Discrete energy inequality). Let (p^, J"^) be an ap- 
proximate solution of the QHD system, with < r < 1. Then, for 
t G [Nt, {N + l)r) we have 

N 

(59) E^{t) <-\Y. l|A(A:r-)|U.(K3) + (1 + t)Eo. 

k=l 

Proof. For all A; > 1, we have 

E^{kT+) - E^{kT-) = J i|A-(A;r+)p - l|A-(fcr-)|Ma; 

J- J{-2r + r'W{kT~)\' 

+ 2(1 - T)A^{kT-) ■ Rk + \Rk\'^dx 

<^ j (-2r + T^)\A^{kT-)\^dx + (1 - T)a\A^{kT- 
1 — r 

+ \Rk? + \Rk?dx 

a 

=^ j (-2r + T^ + a- «r)|A^(A;r-)|2 

+ (^^^) m'dx. 

Here Rk denotes the error term as in espression fl53|l . If we choose 
a = T, it follows 

(60) E-ikT+) - E^ikr-) < -I||A-(fcr-)||i. + ^\\Rk\\l. 

<-^\\A^{kr-)\\l.+r2-'-^MHU 

The inequality ( l59i) follows by summing up all the terms in ( l60i) and 
by the energy conservation in each time strip [kr, {k + 1)t). □ 



Unfortunately the energy estimates are not sufficient to get enough 
compactness to show the convergence of the sequence of the approxi- 
mate solutions. Indeed from the discrete energy inequality, we get only 
the weak convergence of Vip^ in L°°([0, oo); iJ^(M^)), and therefore the 
quadratic terms in (H9l) could exhibit some concentrations phenomena 
in the limit. 
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More precisely, from energy inequality we get the sequence {ip'^} is uni- 
formly bounded in L°°([0, oo); H^{R^)), hence there exists e L°°([0, oo); H^{R^)), 
such that 

ij^^ifj inL°°([0,oo);iJi(M3)). 

Therefore we get 

V^^y^ inL°°([0,oo);/7i(R3)) 
A"^A inL°°([0,oo);L2(M3)) 

r^^A inL°°([0,oo);LL(M3)). 

The need to pass into the limit the quadratic expressions leads us 
to look for a priori estimates in stronger norms. The relationships 
with the Schrodinger equation brings naturally into this search the 
Strichartz-type estimates. The following results are concerned with 
these estimates. However they are not an immediate consequence of 
the Strichartz estimates for the Schrodinger equation since we have to 
take into account the effects of the updating procedure which we im- 
plement at each time step. 

The following Remark [271 will clarify in detail why we stated the Lemma 
[T71 in place of using the exact factorization property provided by the 
Lemma [T5l 

Lemma 25. Let ip^ be the wave function defined by the fractional step 
method (see Section\Bi), and let t G [Nt, {N + l)r). Then we have 

N 

(61) Vr{t) = U{t)V4Jo - ^ J E - (0^A-(A;r-)) 

k=l 

-I / U{t- s)F{s)ds + y^U{t-kT)rl, 

"^0 fc=l 

where as before U (t) is the free Schrodinger group, 

(62) ||0^||loo(k3) < 1, ||r^||L2(K3) < T\\i!o\\Hi{RS) 

and F = V(|V^^|p~V^ + V^ip^). 

Remark 26. Clearly, with a general collision term f{^., J, V^/p) (see 
the Remark{5\) and a non-zero doping profile C{x) one has F defined 
by F = V{\il''^\^~^ip'^ + V^ip'^ + g'^ip'^), where now solves 

-AV^ = |^"|2-C(x) 

and g^ = g{t, x, ^/^, A^, W ^/^) . 
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Remark 27. Now we can see why the updating step in ( ISTl) has been 
defined through the Lemma\r^ Indeed, one could possibly try to use the 
exact unitary factor as in the Lemma lTR namely ip'^ {kT+) := (f)^^"'^^ 
where 0, ^/p are respectively the unitary factor and the amplitude of 
ip'^{kT—) (note that (J) is uniquely determined ^/pdx , and |0| = 1, ^/pdx- 
a.e.) The meaning of this definition is clear if we approximate ip'^ {kT—) 
( as in the proof of the Lemma [7?y with smooth ipn = e*^" ^/p^ and then 
we update ipn withipn '■= e*^^~^^^"y^. Hence ip{kT+) := limn^ooV^n in 
H^{M.^). Moreover in this case we get the exact updating formulas 

\r{kT+)\' = \rikT)\' 

hlm{^Vip^){kT+) = (1 - T)hlm(i^Vip^){kT-), 

without the small error in the expression for the density current, 
however instead of the formula (!6T!) we would find 

Vr{t) =U{t - Nr)a],,U{r) . . . alU{r)V^o 

- i^U(t - NT)al^k^{NT-) + ...- i^Uit -Nt)... U{t)(PIA^{t-) 
n n 

-i / U{t - s)F{s)ds - iU{t - Nt)(tI^^ / U{Nt - s)F{s)ds 

J Nt J{N-1)t 

+ ...- lUit - NT)a],^U{T) ...al [ U{t - s)F{s)ds, 

Jo 

where aj,^ = (0^^)"*, 4>It- being the polar factor of ip'^^kr—). Now, to 
recover an expression similar to fl6ip . we have to calculate the commu- 
tators between the free Schrodinger evolution operator and the multi- 
plication operators by al^. The estimates of these commutators are not 
known for non-smooth al^. 

Proof. Since ip'^ is solution of the Schrodinger-Poisson system in the 
space-time slab [Nt, {N + l)r) x M'^, then we can write 

(63) V^J^it) = U{t- Nt)ViI)^{Nt+) -i ! U{t- s)F{s)ds, 

J Nt 

where F is defined in the statement of the Lemma [251 Now there exists 
a piecewise smooth function 6* at, as specified in the proof of the Lemma 
fTTl such that 

^(Arr+) = e'(i-")^'^v^, 

and furthermore all the estimates of the Lemma [T7| hold, with = 

ip'^i^Nr—), ip = %Ij'^{Nt+) and e = 2~^t\\iPo\\hi{r3)- Therefore, we have 

(64) V^p^(NT+) = V^/j^(Nt-) - ile'^^-^'^'^''A^(NT-) + rL, 

n 
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where ||r]v||i2 < rH-i/^ollHi- By plugging ( IMI) into ( l63l) we deduce 

WVt) = U(t - Nt)V4j^(Nt-) - i^Uit - NrMe'^^-^^'^'^K'iNT-)) 

n 

+ U{t- Nryi, - I [ U{t- s)F{s)ds. 

Let us iterate this formula, repeating the same procedure for V^p'^{NT—) 
then ([61]) holds. □ 



At this point we can use the formula ( 1611) to obtain Strichartz esti- 
mates for Vip'^, simply by applying the Theorem [9] to each of the terms 
in ( pTl) . After some computations, the following result holds. 

Proposition 28 (Strichartz estimates for Vip"^). Let, < T < oo, let 
ip'^ be as in the previous section, then one has 

1 

(65) l|V^^||L^Lr([o,T]xR3) < C(EoM|po||li(R3),T) 

for each admissible pair of exponents {q, r) . 

Remark 29. As we showed in Remark lMj in the general case the non- 
homogeneous term F is slightly different. Anyway by the Strichartz 
estimates the same result of Proposition [13 holds even in the case of 
modified non-homogeneous term F. 

Proof. First of all, let us prove the Proposition [281 fo^' ^ small time 
< Ti < T and let (g, r) be an admissible pair of exponents. We 
choose Ti > later. Let be a positive integer such that Ti < Nr. 
By applying the Theorem [9] to the formula (!6T|) , we get 

ll^''^'^llLt^L5([0,Ti]xM3) ^ ll^(^)Wo|lL9^r([o,Ti]xR3) 

N 

+ lY.\Ht- kr) (e^(^-^)^^A^(A;r-)) |L,,.([o,T,xn 

k=l 

N 



+ Y,\\U{t-kr)rl 
=1 

/ U(t-s)F{s)ds 
Jo 



L?Lj([0,ri]xIR3) 

k=l 

t 



L?LS([0,Ti]xI 



--: A + B + C + D. 



Now we estimate term by term the above expression. 
The estimate of A is straightforward, since 

||f^WV^o|U?LS([0,Ti]xR3) < ||V^/'o||l2(R3). 
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The estimate of B follows from 

N 



29 



(66) I Y: pit - kr) (e^(-^)^^A^(A;r-)) 



xR3) 



k=l 



N 



< 



■5^||A^(A;r-)|U.(M3)<Tii?| 



k=l 



The term C can be estimated in a similar way, namely 



N 



N 



(67) Yl - f^'^yilkLr < Y W^IWlH^^) ^ TMw 



k=l 



k=l 



The last term is a little bit more tricky to estimate. First of all 
we decompose F into three terms, F = Fi + F2 + F^, where Fi = 
VdV'^P'V^), ^2 = VV^ip^ and F3 = V^Vip^. By the Strichartz esti- 
mates (Theorem [H]), we have 



(6J 



[ U{t-s)F{s)ds 
Jo 



L?L-([0,Ti]xR3) 



< 11^1 1 



L^iL,i([0,Ti]x]R3) 



+ m 



Ll^L/{[0,Ti]xR3) 



'L^3i^3([0,Ti]xII 



where {qi,ri) are pairs of admissible exponents. Let us start with the 
first term. 

Lemma 30. There exists a > 0, depending on p, such that 



(69) 



Li LS'([0,Ti]xM3) 



5i([0,Ti]xII 



Proof. First of all let us apply Holder inequality on the left-hand side 
of (!69|) . Then we have 



(70) iWij^fvij^ 



Ll L5'([0,Ti]xR3) 
' r|p-l| 



a||„/,r||p-l 
t 



Now we want 



. — pr = I ( I r— TT ) and — = | ( ^ — — ) , in such a 

way that ||/|Ui{p-i) .n(p-i), l|V/||ipi-2 < H/H^i = HV/H^o. We already 
know ^7 = ^ + \ {\ — jr), then putting together the conditions on 



30 PAOLO ANTONELLI AND PIERANGELO MARCATI 

(g,f), {qj,rj), it follows 

1 1 1 1 1 , ,3/1 
-7 = - + - + - = - + (^^-l)oU 



q' a qi q2 a 2 \6 ri(p— 1) 



^ 2 V 2 rai" ^ 2 1 2 f 



and when 1 < p < 5, 



a = > 0. 

4 



This means that we can always choose pairs (g', f'), (gi, ri), (g2, ^2) sat- 
isfying the previous conditions so that the inequality (!70l) holds, with 
a > 0. 

For instance, if 1 < n < 3, we can choose — = = h, there- 

fore qi = q2 = 00, hence we have ^ = = In the case 

3 < j9 < 5, we take ^ = j~ — ^ then gi = oo,q2 = 2, hence we 
have = I,^ = □ 

r' 6 ' 5 4 



Now let us consider the second term: V'^'Vip'^ , here we choose (^2, '^'2) ~ 
(1, 2), so by the Holder and the Hardy-Littlewood-Sobolev (see [19]) in- 
equalities one has 

\\y^^i^^\\LjLU[0,T,]xR^) < Ti^\\V^\\L^L'i\\^i'^\\L^^L% 



For the third term, we choose (^3, r'^) = (2^^ T+2e^ again by using 
the Hardy-Littlewood-Sobolev and the Holder inequalities, we have 

"Lf=^Li^([0,Ti]xR3) - 111 IIlT^lIII"^ 11"^* 



1 



2 6 
T l+tin 



<T/l|vrni^^||^1U^.2<r,^ll^-|li.,.||vv^1l _ 

Now, we summarize the previous estimates by using (16T|) in the follow- 
ing way 

(71) II VV^1Uo([o,T,]xR3) < II V^o||l^(m3) + Ti4 + II VV'1l5„(j„^^^j^^3) 

1 ^ 

+ II'^o|Il2(k3)||V'?/'^||50([o,Ti]xR3) 

<(l+T)4+Tr||V^1|J„(j„^^^j^^3^+T/||V^o||i2(^3)||V^1Uo([o,T^^ 
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T) > andCi{Eo, ||V^o||l2(m3), 



Lemma 31. There existTi{Eo, \\iPo\\l'2 
0, independent on t, such that 



(72) l|V^1Uo([o,T]xM3)<C^l(i5o, 

for allO<f <Ti{Eo,\\Mlh 



T) > 



T) 



Let us recall that E'^{t) = Eq and ||'?/''^(t)||L2(R3) = ||^/'o||l2(]r3), hence 
we are in the situation in which we can repeat our argument every 
time interval of length Ti, depending always on the same parameters 
Eq, \\iPq\\l2. The consequence of this fact is the following inequality on 
[0,T] 



(73) \\vr\\ 

< Ci(Eo, 



^o||l2,T) 



T 



C{\\Ml^,Eo,T). 



□ 



Proof of the Lemmal3l\ Let us consider the non-trivial case ||^/'o||l2 > 0. 
Assume that X G (0, oo) satisfies 

(74) X < A + fiX + XXP = (f){X), 

with p > 1, y4 > and for all < /i < 1, A > 0. Let X* be such that 

1 

4>'{X^:) = 1, namely X* = (^^^ j ^ \ hence one has 0(X^,) < X^ each 
time the following inequality is satisfied 

(75) 

\pp-i J Xp~i 

Therefore the convexity of (p implies that, if the condition fl75|) holds, 
there exist two roots X±, X+{fi, X, A) > X^, > X^{fi, X, A), to the 
equation (f){X) = X. It then follows either < X < X_, or X > X+. 
In our case /i = Tj'^^^||'?/^o|li2, X = T^, A = (1+T)i?y^, hence we assume 
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Therefore we choose 



(76) Ti := min 



(2||^o||l2^ 



p 



1 p_ 

p-1 — p p-1 



2^(i + r)Eo 



1/2 



4(P-1) 
5-p 



Clearly we cannot have 

'1 - TiIIV^oIIl^ 



1 

p-i 



<x+< WVtp^l 



^"([O.TilxRa)' 



since we get a contradiction as Ti ^ 0, hence 



(77) 



< x_ 



□ 



Corollary 32. Let < T < oo and let ^/]f',^J' he as in previous 
section, then 

1 

(78) ||Vv/p^||l«l-([o,t]xk3) + ||A''||L9^r([o,T]xiR3) < C{E^ , ||po||li(r3),T), 
for each admissible pair of exponents {q, r). 

Unfortunately this is not enough to achieve the convergence of the 
quadratic terms. We need some additional compactness estimates on 
the sequence {Vip^} in order to apply Theorem [T3l In particular we 
need some tightness and regularity properties on the sequence {'Vip'^}, 
therefore we apply some results concerning local smoothing due to Vega 
and Constantin, Saut [9J. 

Proposition 33 (Local smoothing for V^^"^). Let < T < oo and let 

ip'^ he defined as in the previous section. Then one has 

(79) \\vr\\L.(^^,^T];H]i:m) - ^(^0'^' ii/^oii^o- 

Proof. Using the Strichartz estimates obtained above, we can apply the 
Theorems [TTl [12] about local smoothing. Indeed, by using again the 



FINITE ENERGY WEAK SOLUTIONS TO QHD 33 
formula ( 16T1) it follows 

ll^^"llL2([0,T];<e(M3)) <l|V^o|U2(R3) 

N 

k=l 

N 

+ r^||V^„JU2 

k=l 

N 



fc=i 



L2 



+ l|-^l|Li([0,r];L2 



The first three terms are clearly estimated by a constant C{Eo,T) 
depending only on the initial energy and on time. The fourth term is 
0{t). The last term will be estimated using the previous Strichartz 
estimates. As before, we split F into three parts F = Fi + F2 + F^, 
then one has 

\\W\^'^^V\\lIlI([0,T]xR3) 

<T^\\\^Ij''\p~'^\\ 4 \\^^''\\l°-lU[o,t]xr3) 

<^^ll^1li7L^([0,T]xR3)l|V^1U-L2([o,T]xR3), 

while 

and now the remaining calculations are similar to those already done 
for the Strichartz estimates. Regarding the term V'^'Vip'^ we already 
estimated its LjL^ norm. □ 



Since H^^^ is compactly embedded in Lf^^, we can apply the Theorem 
[T3] due to Rakotoson, Temam [H] . 

Proposition 34. The sequence {'Vip'^} is strongly convergent in 
L\[0,T]-LUR^)), namely 

(80) W := s - lim V^^" m ,T\, L]^J^)) . 

In 'particular, onehas^J ^ \J ^ and A m L^([0, T]; Lf^^(M^)). 
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Proof. The previous Proposition [331 implies that the sequence {VV^^}t->o, 
is uniformly bounded in L^([0, T]; if//^^(M^)) and then Vip'^ Vip in 

that space. Now H^^^ is compactly embedded in Lf^^, since 'Vijj'^{t) 
Vip{t) for almost every t > and 

lini^ _sup / ||V^^(t)||^2 dt = 0, 



\E\^0,EC[0,T] r>0 Je 



loc 



since Vi/j'^ G L°°([0, T]; L^(M^)). Hence we can apply the Theorem [T31 
by Rakotoson and Temam [H] and we get flHU]) . □ 



Proposition 35. (p, J) is a weak solution to the Cauchy problem ([T]), 
0. 

Proof. It follows directly by combining the Theorem [251 in the section [S] 
and the Proposition [311 As for the collisionless QHD system we should 
note that the generalized irrotationality condition holds by the defini- 
tion of the current density and Corollary [161 D 
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